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Abstract 



Explicit expressions for the expectation values and the variances of some 
observables, which are bilinear quantities in the quantum fields on a D- 
dimensional manifold, are derived making use of zeta function regularization. 
It is found that the variance, related to the second functional variation of the 
effective action, requires a further regularization and that the relative reg- 
ularized variance turns out to be j^, where is the number of the fields, 
thus being independent on the dimension D. Some illustrating examples are 
worked through. The issue of the stress tensor is also briefly addressed. 
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I. INTRODUCTION 



Vacuum fluctuations play an important role in many physical processes. The Casimir 
effect is one of its most interesting physical manifestations and it has been experimentally 
verified. It is also well known that the Casimir effect is related to the presence of a non 
vanishing vacuum energy (see for example This fact mainly occurs when one is dealing 

with non trivial space-times, where the spatial sections are topologically non trivial spaces 
or manifolds with boundaries. 

Another interesting issue where quantum vacuum effects are present is the Quantum 
Field Theory in curved space-time Recall that within the semi-classic approach to 

quantum gravity, the basic equation reads 

G^^ + g^^X = SttG < T^^ > , (1) 

where G^jy is the Einstein tensor, A the cosmological constant and T^i, is the vacuum expec- 
tation value of the matter stress tensor (we use h = c = 1). As a consequence, fluctuations 
of the stress tensor can induce fluctuations of the classical gravitational field and, in order to 
justify the semiclassical approximation, it appears very important to have reliable a priori 
estimates of these fluctuations. 

Fluctuations of the stress tensor were studied in making use of canonical methods. 
Fluctuations of Casimir forces were investigated in Alternatively, other authors |^ 

have investigated the same problem by making use of zeta function regularization. 

With regard to this issue, it is well known that the notion of effective action (or effective 
potential) plays an important role as a powerful tool in relativistic quantum field theory. 
This quantity, however, is ill defined, since, within the Euclidean formulation and in the 
one-loop approximation, the one-loop effective action contains functional determinants of 
elliptic operators, which have to be regularized. Zeta-function regularization [1T0|-|1^, (see 
also |T^ for the generalization to elliptic pseudo-differential operators and p5|-p!7| for physical 



applications) was introduced by a number of authors as a convenient tool in order to deal 
with the evaluation of functional determinants. It permits to give a meaning in the sense 
of analytic continuation — a mathematically very precise procedure — to quantities that are 
formally divergent. 

In this paper, we would like to revisit the zeta function regularization approach for the 
evaluation of expectation values < O > and their quantum fluctuations. It is our opinion 
that zeta function regularization is a very powerful tool at our disposal, as compared with 
other methods and that this issue deserves a careful investigation. 

We will mainly consider two quantities: O = cfP' and the stress tensor trace O = Tj^ 
and their corresponding variances A^O =< > — < O >^. Within our formalism, it is 
convenient to introduce the relative variance ^ 

A^0 = ^^. (2) 

A different relative variance, though directly related to the previous one, has been introduced 
by Kuo and Ford 0. It reads 
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A2/9 A2/9 

Considering now the operator O = Tj^ , we observe that the trace of the Einstein equations 
reads 

G''^ + DA = 87iG <TI^ > . (4) 

Thus, fluctuations of the stress tensor trace induce classical fluctuations of the scalar curva- 
ture. Furthermore, this trace fluctuation may be used to have an estimate of the fluctuations 
of the whole stress tensor, since for conformally invariant quantum flelds in homogeneous 
space-times, one has 

< T,. >= ^ < > . (5) 
The issue concerning the validity of semiclassical gravity has been discussed in []6|,|T8 



We also recall that the flrst variation of the effective action is related to the vacuum 
expectation value of physical quantities, while the second variation of the effective action is 
associated with the variance. Within the zeta function regularization procedure, one has to 
deal with traces of complex powers of elliptic operators. The flrst variation of the effective 
action is well deflned by the use of zeta function regularization, while the second variation 
is intrinsically ill deflned, unless one makes use of suitable variations with disjoint supports. 
In the coincidence limit, the physically interesting case, one has to make use of a further 
regularization 0. 

Our main result is the following: modulo regularization problems, the relative variance 
turns out to be A^O = j^, (thus A^O = j^), where N is the number of scalar flelds in 
some multiplet. This result seems to be general, that is independent on the quantity one is 
dealing with, for example the stress tensor trace. Our results are compatible and should be 
compared with the ones recently obtained, for N = 1, regarding the vacuum energy density 



fluctuations via smeared flelds and point separation [|T9| , which give relative variances of the 
order of unity albeit dimensionally dependent. 

On the other hand, coming back to the case of neutral scalar non-interacting flelds, we 
recover a well known criterion for neglecting the quantum gravity fluctuations in the large 
limit 13 . 

The content of the paper is as follows. In Sect. II, zeta function regularization and heat- 
kernel techniques are briefly summarized. In Sect. Ill, the flrst variation of the effective 
action computed is shown to be related to the vacuum expectation values of observables. 
In Sect. IV, the second variation is shown to be associated with the variance and the flnal 
result is presented with the help of a further analytic and ad hoc regularization. In Sect. 
V, some examples are presented. The paper ends with some concluding remarks and an 
Appendix, where the flrst and second variations of the trace of the complex power of an 
elliptic operator are explicitly computed. 



II. ZETA FUNCTION REGULARIZATION OF THE EFFECTIVE ACTION 

In this section, we will summarize some basic aspects of the heat-kernel and zeta function 
regularization methods. For the sake of simplicity we shall here restrict ourselves to scalar 
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fields, but the method is also valid in more general situations. In the case of a neutral scalar 
field, the one-loop Euclidean partition function, reads [12| 



= -InZ = -lndet4 , (6) 
2 /i^ 

where /i^ is a renormalization parameter. 

To begin with, recall the definition of the zeta function related to an elliptic operator A 

1 

mas\A) = -— f-'Tr e~'^dt, (7) 
i [s) Jo 

which is valid for large values of 3?s. 

For small t, and for second order elliptic differential operators, it will be assumed that 
the heat-trace has the following asymptotics 

oo 

TTe~'^ = j:K4-^, (8) 

r=0 

where are the integrated heat-kernel coefficients. In principle they may be computed 
(see, for example [0). We will also assume the validity of the local heat-kernel expansion, 
namely 

oo 

e-'''{x) = J2Kr{x)t'^. (9) 

r=0 

If the asymptotic expansions and (|^) hold, a standard argument leads to the following 
meromorphic extension of the zeta function and its local counterpart {J{s) being an analityc 
function) 

ns)as\A) = ^^j^ + j{s). (10) 

r=0 + 2 



T{s)as\A){x) = E -^TED + As, x) . (11) 

r=0 S + 2 

As a result, at s = 0, the global and local zeta functions are regular and their derivatives 
exist. We also have C,{Q\A){x) = Kj:,{x). As is well known, Kj:,{x) = in any odd dimen- 
sional manifold without boundary. Under hypothesis above, zeta function regularization 
corresponds to doing the following: 

,„detA = _c-(0|^, = -,tolT.(,-^-.,. (12) 

However we would like just to recall that, in some situations, it might be necessary to 
generalize the above definition in the form 

Indet 4 = -lim i-^Tr(./.2M-) . (13) 
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When Tr(yl~*) is regular at s = 0, the two definitions (|T^ and ([131) coincide, but, in some 
cases, C("S|^) has a simple pole at s = and so Eq. (12) is no more valid. This fact may 
well be present at the level of the effective action (see for example [p^-p^). 

For the sake of simplicity, in the sequel we will assume the validity of Eqs. ([T0|) and (p|). 
The latter gives 

mA){x) = Kd{x) . (14) 

For future use, we also observe that when K£,_^2{x) 7^ 0, the local zeta function has a 
simple pole at s = 1. In fact we may write 

C(l + s\A){x) = ^^^^ ^^^^ + G(l + s,x), (15) 

where G{l,x) is a regular function. It is given by 

G{l,x) = PPai\A){x)-jKn-2{x) , (16) 
7 being the Euler's constant and 

PPC(l|A)(x) = Ihn (^C(l + s\A){x) - ^^^iM^ , (17) 

the finite part of C('5|^) at s = 1. When we consider an odd dimensional manifold without 
a boundary, this singularity is absent. 



III. THE FIRST VARIATION OF THE EFFECTIVE ACTION: VACUUM 

EXPECTATION VALUES 

In this section, we will evaluate vacuum expectation values < O > of some specific 
quantities O, such as the stress tensor trace or conformal anomaly and the square of the 
field (fiuctuation). These quantities involve the product of two quantum fields at the same 
point, X, and are therefore ill defined. They require a regularization. We shall consider a 
multiplet of N scalar fields denoted by (p in an external field, described through a classical 
action given by 

S=^Jdx <j)L(t) , (18) 

where L is a suitable (matrix valued) differential operator defined on a D-dimensional smooth 
manifold. 

To begin with, let us recall the formal trick that allows one to get the vacuum expectation 
value of the bilinear O = (pKcf) within a path integral approach |^|. We will consider 



two cases: K = J, the identity matrix, in the case of the field fiuctuation O = (f)"^, and 
K = Ci + C2L, in the case of stress tensor trace O = T^^. Here, Ci and C2 are constants and 
moreover, in the conformally invariant case, Ci = 0. If we denote by a{x) a suitable classical 
source, we may consider the Euclidean partition function 
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Zia) 



-^Jdx^Lia}<f> 



det 



L{a) 
1^ 



-N/2 



(19) 



Here we are assuming that, in the massive case, the multiplet has the same common mass. 
In this way, there is no multiphcative anomaly (see, for example, and L{a) = L + aK 
may be regarded simple differential operator. 
A formal functional derivation leads to 



< 0{x) >-- 



6lnZ(a) 



6a 



(20) 



a=0 



Making use of (|TU[) and zeta function regularization (see Eq. ( [55| ) in the Appendix), we may 
give a meaning to the above formal expression, namely 



N 



d 



In the case O 



S In Z(a) = lim — 

6L{a) = 6a. Thus 



< (f)^(x) >= A^lim-^ 
^ ^ s^o ds 



/i^^sTr (L-'-\a)6L{a 



s^i^X{s + l\L){x) 



(21) 



(22) 



In the equations above, C,{s\L){x) is the local zeta function. As a result, making use of the 
meromorphic expansion (|15]), one gets (see |2^^^9|) 



< 02(3,) A^PPC(l|L)(x) + NKd-2{x) In/i^ . 

When D is odd and the manifold is without boundary, we simply have 

<cP\x) >= Nai\L){x) . 

In the other case, namely when one is dealing with the stress tensor trace, 6L{a) 
6a{ci + C2L). As a consequence. 



(23) 



(24) 



< Tl'ix) >= Nlim^- 



sfi'^cCi-s + l\L){x) + C2Cis\L){x) 



and, result, 



< Tl^{x) >= ci < 0^(x) > +C2NKd{x) . 



(25) 



(26) 



In the conformally coupled case Ci = and one has the usual conformal anomaly, due only 
to quantum effects. 
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IV. THE SECOND VARIATION OF THE EFFECTIVE ACTION: THE 

VARIANCE 



We have seen that the first variation of the effective action is associated with the vacuum 
expectation value < O > of bihnear quantities in quantum fields. Let us show now that the 
second variation of the effective action is related to the variance A^O =< > — < O >^. 

To begin with, the second variation of the partition function ( [ToD gives 



< 0{x)0{y) >= 4- 



1 6^Z(a) 



Z{a) 6a{x)6a{y) 
On the other hand, we have the identity 

6HnZ(a) 1 6^Z(a) 



(27) 



a=0 



1 SZ{a) SZ{a) 



6a{x)5a{y) Z{a) 5a{x)5a{y) Z'^{a) 5a{x) Sa{y) 
As a consequence, the variance is given by 

6^ In Z(a) 



(28) 



A'0{x,y) =< 0{x)0{y) > - < 0{x) >< 0{y) >= 4 

Now, it appears convenient to introduce the relative variance 

A^Oix,y) 



6a{x)6a{y) 



(29) 



a=0 



^:Oix,y) 



< 0{x) >< 0{y) > 



(30) 



The coincidence case, x = y, is particularly interesting from the physical point of view. It 
is formally given by 



AlO{x) = AlO{x,y 



\x=y 



6^ In Z{a) f 6lnZ{a) 
5'^a(x) \ 6a(x) 



(31) 



Let us evaluate the second variation of lnZ{a). In general, in the coincidence limit 
one gets an ill defined quantity and a further regularization is required, as explained in the 
Appendix. Making use of (P^)? (0) (PD the Appendix, one has 



62S1 In Z{a){e) 



2e 



lim — 



2s 



T s) 



2r2(l + e) s^ods 

poo roo , 

duu' j rfww=(M + t;)^Tr (e-"%Ae-^^5iA 
For e > and sufficiently large, the integrand is regular at s = 0, and we have 



(32) 



52^1 \nZ{a){e) 



/* 00 /* 00 

duu' dvv'Ti [e~''^62Ae-''^6iA 



212(1 + e 

.2s 



(33) 
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This is our general formula for the second variation of the functional determinant. It is in 



agreement with a result obtained in |T3 



Let us consider the specific variation related to our observables. First, for O 
SL{a) = Sa, and in the coincidence limit we have 



< A'(j)' > (e) = 2Npi'X\l + e\L){x). (34) 

For X 7^ y, (,{1 + e\L){x,y) is regular at e = 0, but in the coincidence limit y ^ x, C{l\L){x) 
is well defined only for odd D and boundary free manifolds. 

In the conformally coupled case and for O = T^, one has 5L{a) = C2SaL. As a result 

< A^T; > (s) = 2Ncie{e\L){x) = 2iVc^C'(0|L)(x) + 0(5) (35) 

and we can remove the regularization parameter because (^{s\L){x) is regular at the origin. 

Some remarks are in order. If we limit ourselves to the odd dimensional case, in a 
manifold without boundary, the analytic continuation works in a simple manner for both 
the quantity O = and = 0^ and no scale dependence appears in the final expressions. 
As a consequence, it turns out that the relative variance, 

A^O = , (36) 

is always equal to 

A^O = I , (37) 

and 

AlO = . (38) 

N + 2 ^ ' 

The expressions above are valid for a multiplet of scalar fields. In the case of other 
quantum fields, the formula for the relative variance could be a little bit different (see for 
example ref . for the treatment of gauge fields) . 



In the even dimensional case, or/and in the presence of boundary, the situation is more 
complicated and a further renormalization seems unavoidable. 



V. EXAMPLES 
A. D— dimensional torus 

In this Section we will consider, as first example, the D— dimensional torus. This is a 
symmetric flat manifold with finite volume and the local zeta function is simply the ratio of 
the global zeta function and the torus volume. Thus, we may limit ourselves to the discussion 
involving the zeta function. This is a general conclusion valid for every symmetric space. 

The zeta function for the case of a massive (A^ = 1) scalar field is given by 
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Cd(s|L) 



(39) 



with 30 - 33 



S/2-D/4 



r{s-D/2) 2"/^+^/^+V^(mi?)-^+^/^ 



Kd/2-s h?''^TxmRV^ 



Vis) 



(40) 



Here R is the radius of the torus and m the mass of the field. This expansion is exponentially 
convergent. It is to be seen not just as a big mass expansion (the convergence is then 
extremely fast), for it is valid in a very wide range of values of mR: 1 ~ mR < oo. 

In the case of a massless field, the convenient expression to be used is quite different 
(this is explained in detail in [Q): 



2l+s D-l 



■^■/2r(s-j72)C/j(2s-j) 



n 



i/2- 



n=l 



s/2-j/i 



(41) 



On the other hand, for the sake of completeness, it is interesting to have a perturbative 
expression for the case when the mass of the field is very small but different from zero. This 
is obtained by means of binomial expansion in the equation defining the zeta function. The 
result is: 



E( 

k=0 



[s + k). 



(42) 



This expression, combined with the preceding one, Eq. (|41|) , yields the desired low mass 
expansion. Such expansions for the zeta function are not to be found in the literature. In 
fact, explicit expressions of the type 



yj) have appeared in the seminal paper for the 
first time. They are convenient, on the one hand, because they exhibit the pole structure of 
the zeta function explicitly. On the other, they are useful from the computational point of 
view, because they consist of a term including the main contribution, together with a series 
that converges extremely fast: only a few first terms need to be computed in order to obtain 
results as accurate as desired. 

From these expressions, the first and second variations of the effective action are imme- 
diate to compute. Essentially, what we get are expressions of the following kind: (i) in the 
massless case {p is a natural number, a small one for any of the operators here considered) 



Zd{p) 



[Zd{p-1W 2P-^T{p) 



X 



(43) 



{Ef=V k^/'r(p - 1 - j/2)Cr{2p - 2 - j) + ^^'SoAp - 1)]} 



2 ' 
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where the Sd j{p) are fast convergent series providing only corrections to the main terms, 
fii) 



Zd{p] {mRf 



[ZD{p-l-{mRY)Y 
Zd{p) 



(44) 



[Zd{p-iW 



l + 2p 



Zpjp+l) 

\Zoip-iW [ZDipW 



Zd{p) 



[mRf + 



in the case of a field of very small mass (after doing a small-mass expansion as described 
above), and (iii) 



Zd{p- {mRf 



(m 



RfP~^T{p 



[Zd{p - 1; (mi?)2)]2 (27r)^/2r(p) 

T{p -D/2) + 2P/2-^/4+i7rP-^/2^mi?)P-^/2 + Sz)(p) 



(45) 



[T{p -l-D/2) + 2P/2-i5/4+i/27rP-i-^/2(mi?)p-i-^/2 + j:^(p - i)] 



2 ' 



in the massive case, where again Tj^ip) is a very fast convergent series. These expressions 
simplify very much when poles of the gamma function appear (for even dimension D). They 
are quite easy to deal with. As advanced before, no singularity appears for D odd, and the 
whole expression remains, in that case (the series being well approximated by a couple of 
terms) . 



B. The Casimir slab 

As a second example, we will revisit the computation of the local zeta function related 
to the Casimir slab, namely a massless scalar quantum field confined in the x direction 
between infinite parallel planar Dirichlet boundaries, located at a; = and x = a. In this 
case, the local zeta function is not trivial, due to the presence of the planar boundaries (see. 



for example P,|3|, p!5| , p!6| ) . 

To start with, recall the form of the local diagonal heat-kernel, which depends only on 
the confining coordinate x, i.e. 

. _.r. 2 ~ . ^nnx n^Ti'^t. . 
< x\e \x >= ^ 2^ sm^ exp —) . (46) 



a{ATTt)' 



2 



n=l 



a 



Here L = — is the Dirichlet Laplacian in the slab. Making use of the trigonometric 
duplication formula for the sine and Poisson-Jacobi re-summation formula, we may rewrite 
the above expression in the form 

I -tLi ^ f ^ I I (na + x)2 \ 
<x\e V>=—-n\ E exp( — )- E exp(-^ • (^^) 



(47rt) 



-oo 



All the Seeley-DeWitt K^ix) coefficients vanish, but the first one Kq[x) = 1. As a conse- 
quence, we may anticipate that < T^[x) > is vanishing. Making use of the Mellin transform 
and the above expression for the heat-kernel, and analytically regularizing the integral 
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poo 

/ dte = 0, 
Jo 



(48) 



the analytic continuation for the local zeta function may be obtained. We present the result 
in a simple and symmetric form (for another equivalent form corresponding to D = 4, see 

m 



C{sL){x) = — X 

(47r)fr(s) 

[2Cn{D - 2s) - Ch{D - 2s, x/a) - Ch{D - 2s, 1 - x/a)] . 



(49) 



In the above expression, C,r{z) ^md CH{z,q) are the Riemann and Hurwitz zeta functions 
respectively. This result can also be obtained by making use of the re-summation techniques 
explained in |35| , p!5[ . 



Some comments are in order. The local expression we have obtained is already in the 
renormalized form. Furthermore, the related trace involves a (infinite, but trivial) volume 
Vt in the transverse directions and an integration over x. Since 



/ dq CHiz,q) = 0, 
Jo 



the resulting zeta function reads 

C(^l^) 



2VtT(^ 



s a 



2S-D+1 



(An 



<r{D - 2s) 



(50) 



(51) 



and this is exactly the zeta function associated with the Casimir slab configuration. 

As far as the application to the evaluation of the vacuum expectation value is concerned, 
first, let us consider D > 2. The zeta function is regular at s = 1 and the result is 



< 0^(x) >= 



r(f 



Da 



2-D 



(An) 



X 



(52) 



[2CRiD -2)-Ch{D- 2, x/a) - Ch{D - 2, 1 - x/a)] 



The expression is finite everywhere and gives a vanishing result at x = and x = a. Due to 
the simple geometry of the planar boundaries, the boundary divergences are not present. 

Things are different in D = 2. In this case, the zeta function has a pole at s = 1 and 
the dependence on the scale appears. For the sake of completeness, we give the result. It 
reads 



< 



[x] >-- 



1 

2^ 



7 + In h In sm — 

vr a 



(53) 



In this case, the boundary divergences are present. 

A few other exact analytic results obtained via zeta function regularization can be found 

in 
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VI. CONCLUDING REMARKS 



In this paper, we have revisited the use of zeta function regularization approach to the 
evaluation of expectation values of physical quantities and their related quantum variances. 
The former are associated with the evaluation of the first variation of the effective action, 
while the latter are related to the second variation. We have shown that the first variation can 
be regularized by use of zeta function techniques, and explicit expressions for the vacuum 
expectation values have been exhibited. For issues concerning the second variation, in 
general, zeta function regularization works well only when one is dealing with off diagonal 
terms, since the coincident limit is highly singular. A further analytic regularization has to 
be introduced to treat the coincidence limit. 

We may summarize our results as follows. For a multiplet of N scalar fields and O = Tl^, 
in the conformally coupled case, analytic regularization gives a relative variance exactly 
equal to j^, independently on the dimension of the boundary free manifold. For O = (f)^, 
we have obtained again a relative variance exactly equal to ^, but now limited to the odd 
dimensional case without boundary. These restrictions can be removed in some particular 
situations, like the case of the Casimir slab, where one is dealing with a flat manifold with flat 
boundaries. The even dimensional case seems to require however further renormalization. 

The formalism can be directly applied to the expectation value of the stress tensor (see 
0]), but only after much more effort and work. Again, the first variation of the effective 
action with respect to the metric tensor gives < T^^, >. With regard to this issue, for an 
evaluation using local zeta function regularization see pE|J5B[ |. 



The second variation is related to the variance, and for the off diagonal case, there is 
again factorization. In the special case of homogeneous space-times and for a conformally 
coupled scalar field, one recovers the O = case. 

VII. APPENDIX 

In order to compute the variation of the trace of an elliptic invertible operator A, one 
has to take into account the fact that the variation (deformation) of A does not commute 
with A. From A~^A = I, we have 

6A-^ = -A-^6AA-^ . (54) 

For the calculation of complex powers of A we can use the Cauchy-Dunford representation 

A-' = — - f dzz^'{A-z)-\ (55) 



27ri Jc 

in which C is a suitable contour on the complex z— plane. As a consequence, we have 

5A-' = ^ [ dzz-\A - z)-^5A{A - z)-^ . (56) 
zvrz Jc 

Making use of the two representations (for Re z > and Re s > 0): 

poo 

T{s){A-z)-'= / dtf-^e'^e-'^ , (57) 
Jo 
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and 

1 1 



dwi-wy'e-'" , (58) 



r(2;) 2iTi Jc 
one gets 

6TtA^' = -sTt{A-'-^6A) . (59) 

If the operator A is self-adjoint, then there exist eigenvalues and eigenvectors, A„ and 
such that v4\E'„ = A„\E'„, and we have 

6TtA-' = — I dzz-'Ti{{A - z)-^5A) 
2Txi Jc 

= -sY.K"^ < i'n.SA'l,, >= -sTiiA-'-^SA) . (60) 

n 

For the second variation, making use again of ( p7| ) and (|55D, one gets 

62SiTtA^' = -sTt{62A-'-^6iA) 

= -s—f dzz-'-^TT((A-z)-^62A(A-z)-^6iA). (61) 
2m Jc 

This expression is valid only if Supp^i has void intersection with Supp52- Since we are 
interested in the coincidence limit y = x, we have to introduce an additional analytic 
regularization. The simplest one is to replace {A — z)~^ with [A — 2;)"^"^ and employ 
the representation in terms of the Mellin transform: 

T{l + e){A- z)-^-' = rfwM^e^^e-"^, (62) 

Jo 

in Eq. ( |63D and then use again ([58|). As a result, we arrive to the formula 

626iTTA-%e) = -— --— / duu' / dvv^u + vYTi (e-^^^aAe-^^^i^) . (63) 

^"'(l + e) T{s) Jo Jo ^ ' 

This representation for the second variation of the trace of a complex power of an elliptic 
operator A is in agreement with the one obtained in making use of a different method 
based on Schwinger's perturbative expansion. 



VIII. ACKNOWLEDGMENTS 



We are grateful to Valter Moretti for a useful discussion. This investigation have been 
supported by the cooperation agreement INFN (Italy)-DGICYT (Spain). EE has been 
financed also by DGI/SGPI (Spain), project BFM2000-0810, and by CIRIT (Generalitat de 
Catalunya), contract 1999SGR-00257. 



13 



REFERENCES 



N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Space (Cambridge Univer- 
sity Press, Cambridge, 1982). 

A. A. Crib, S.G. Mamayev and V.M. Mostepanenko, Vacuum quantum Effects in Strong 
Fields (Friedmann Laboratory Publishing, St. Petersburg, 1994). 

S.A. Fulling, Aspects of Quantum Field Theory in Curved Space-time (Cambridge Uni- 
versity Press, Cambridge, 1989). 

R.M. Wald, Quantum Fields in Curved Spacetime and Black Hole Thermodynamics 
(Chicago University Press, Chicago, 1994). 

L.H. Ford, Phys. Rev. 47, 1 (1993); Phys. Rev. 51, 1962 (1995); 
C.I. Kuo and L.H. Ford, Phys. Rev. 47, 4510 (1993). 

G. Barton, J. Phys. A: Math. Gen. 24, 991 (1991); J. Phys. A: Math. Gen. 24, 5533 

(1991) . 

C. Eberlein, J. Phys. A: Math. Gen. 25, 3015 (1992); J. Phys. A: Math. Gen. 25, 3038 

(1992) . 

N. G. Phillips and B.L. Hu, Phys. Rev. 55, 6123 (1997). 

D. B. Ray and I.M. Singer, Advances in Math. 7, 145 (1971). 
J.S. Dowker and R. Critchley, Phys. Rev. D13, 3224 (1976). 
S.W. Hawking, Commun. Math. Phys. 55, 133 (1977). 

M. Kontsevich and S. Vishik. Functional Analysis on the Eve of the 21st Century. 
volume 1, 173-197, (1995). 

Kirill A. Kazakow. Quantum fluctuations of effective fields and the correspondence prin- 
ciple, |hep-th/02072T8| (2002). 

E. Elizalde, S. D. Odintsov, A. Romeo, A. A. Bytsenko and S. Zerbini, Zeta Regulariza- 
tion Techniques with Applications (World Scientific, Singapore, 1994). 

E. Elizalde, Ten Physical Applications of Spectral Zeta Functions (Springer Verlag, 
Berlin, 1995). 

A. A. Bytsenko, G. Cognola, L. Vanzo and S. Zerbini. Phys. Rep. 266, 1 (1996). 

B. L. Hu and N. G. Phillips, Int. J. Theor. Phys. 39, 1817 (2000). 
N.G. Phillips and B.L. Hu, Phys. Rev. 62, 084017 (2000). 

E. Toumboulis, Phys. Letters B70, 361 (1977); C.G. Callan, S.B. Giddings, J.A. Harvey, 
and A. Strominger, Phys. Rev. 45, 1005 (1992); D. Binosi, V. Moretti, L. Vanzo and S. 
Zerbini, Phys. Rev. 59, 104017 (1999). 

T.P. Branson, P.B. Gilkey, K. Kirsten and D.V. Vassilevich, Nucl. Phys. B563, 603 
(1999). 

A. A. Bytsenko, G. Cognola, and S. Zerbini, Class. Quant. Grav. 14, 615 (1997). 
A.A. Bytsenko, A.E. Goncalves, and S. Zerbini, Mod. Phys. Lett. A16, 1479 (2001). 
A.A. Bytsenko, E. Elizalde, and S. Zerbini, Phys. Rev. D64, 105024 (2001). 
G. Cognola and S. Zerbini, Mod. Phys. Letters 3, 599 (1988). 

E. Elizalde, L. Vanzo and S. Zerbini, Comm. Math. Phys. 194, 613 (1998); E. Elizalde, 
G. Cognola and S. Zerbini, Nucl. Phys. B 532, 407 (1998); E. Elizalde, A. Filippi, L. 
Vanzo and S. Zerbini, Phys. Rev. D 57, 7430 (1998); 
D. leUici and V. Moretti, Phys. Lett. B425, 33 (1998). 
V. Moretti, Commun. Math. Phys. 201, 327 (1999). 



14 



[29] D. Binosi and S. Zerbini, J. Math. Phys. 40, 5106 (1999). 

[30] E. Elizalde, J. Phys. A34, 3025 (2001). 

[31] E. Ehzalde, J. Comput. AppL Math. 118, 125 (2000). 

[32] E. Ehzalde, Commun. Math. Phys. 198, 83 (1998). 

[33] M. Bordag, E. Ehzalde and K. Kirsten, J. Math. Phys. 37, 895 (1996). 

[34] I.M. Gclfand and G.E. Chilov, Les Distributions, tome 1 (Dunod, Paris, 

[35] A. Actor, Fortschr. Phys. 35, 793 (1987); Fortschr. Phys. 43, 141 (1989). 

[36] V. Moretti, Journ. Math. Phys. 40, 33843 (1999). 



15 



